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Abstract 

Recent experiments have demonstrated that the performances of organic FETs strongly depend on the dielectric 
properties of the gate insulator. In particular, it has been shown that the temperature dependence of the mobility 
evolves from a metallic-like to an insulating behavior upon increasing the dielectric constant of the gate material. 
This phenomenon can be explained in terms of the formation of small polarons, due to the polar interaction 
of the charge carriers with the phonons at the organic/dielectric interface. Building on this model, the possible 
consequences of the Coulomb repulsion between the carriers at high concentrations are analyzed. 
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1. Introduction 

In organic field effect transistors (FETs), charges 
move in a conducting layer located at the inter- 
face between an organic crystal and a gate in- 
sulator. Commonly used gate insulators such as 
Si02 or AI2O3 are polar materials. When an elec- 
tron moves in the vicinity of such materials, it in- 
duces a long-range polarization that modifies its 
physical properties, forming a polaron. This mecha- 
nism is well known in inorganic semiconductor het- 
erostructures such as GaAs/AlAs.[l] There, how- 
ever, due to the weak polarizabilities and to the ex- 
tremely low band masses, the polar coupling only 
causes a small renormalization of the electronic 
properties. 



In organic FETs, there are two important differ- 
ences that can eventually lead to much more dra- 
matic effects. The first is that organic crystals have 
very narrow bandwidths, due to the weak Van der 
Waals bonding between the molecules. The second 
is that one can use gate insulators where the static 
(e s ) and high frequency (too) dielectric constants 
have a large difference in magnitude, i.e. that are 
much more polarizable than GaAs. 

Indeed, recent experiments performed on 
rubrene single crystal FETs have shown that the 
charge carrier mobility at room temperature de- 
creases upon increasing the dielectric constant of 
the gate material. [2] At the same time, the tem- 
perature dependence of the mobility evolves from 
a metallic-like to an insulating-like behavior, [3] 
which gives strong support to the relevance of 
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polaron formation. 

In this work we present a theoretical framework 
that consistently explains the observed experimen- 
tal behavior. The model which describes the inter- 
action of the charge carriers with the ionic polariza- 
tion at the interface, as well as the basic results con- 
cerning the formation of small polarons, are pre- 
sented in Section 2. The temperature dependence 
of the polaronic hopping mobility is calculated in 
Section 3, where we briefly discuss the experimen- 
tal results of Ref. [3] . The possible relevance of the 
long-range Coulomb repulsion between the carri- 
ers at sufficiently high concentrations is discussed 
theoretically in Section 4. 



2. Polaron formation 

The interaction between the carriers and the po- 
lar interface is described by the following Frohlich 
Hamiltonian [1,4,5] 

H = e fe c+c fe + frw a ^ a t a i + 

k q 

+ E M 9 c fe +9 c feK + (i) 

kq 

where c^T, Cfe, a+, a q are respectively the creation 
and annihilation operators for carriers, whose dis- 
persion is 6^, and for optical phonons of frequency 
io s . The electron-phonon matrix element is defined 
as 

2 _ 2-Khuj s e 2 (3 _ 2qz 



aS q 
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where a is the lattice spacing in the organic crys- 
tal, S is the total surface and z is the mean dis- 
tance between the carriers and the interface. The 
strength of the electron-phonon coupling is con- 
trolled by the dimensionless parameter 







(3) 



(e s + K)(eoo + k) 

which is a combination of the dielectric constant 
of the organic crystal (k) and those of the gate 
material. 

The above matrix element was derived in the 
long- wavelength limit. A cut-off at short distances 



(at q ~ it /a) is needed in principle to account for 
the discrete nature of the molecular crystal. [7,8] 
However, the finite distance z already acts as a 
short-distance cut-off: M q is exponentially reduced 
at the Brillouin zone boundaries for z > a, in which 
case we can use Eq. (2) for all practical purposes. 
A prescription to deal with discrete lattice effects 
for z < a is proposed in the Appendix. 

Due the interaction term in Eq. (1), the elec- 
trons get "dressed" by the ionic polarization of 
the gate material. If the interaction is sufficiently 
strong, the carriers become self-trapped on indi- 
vidual (or few) organic molecules, forming small 
polarons, and hopping-like transport can set in. 

The polaron energy in the strong coupling 
regime is given by Ep = Y] q M q /u> s . \6] Integrat- 
ing this expression over the whole reciprocal space 
leads to the simple result 



Ep = 2^ 



(4) 



The polaron energy turns out to be independent on 
the phonon frequency ui s , and therefore depends 
on the particular gate dielectric only through the 
parameter (3. It should be noted that the exponen- 
tial decay of the matrix element M q does not im- 
ply an analogous decay of the polaron energy with 
distance, but rather yields a 1/z behaviour as can 
be seen from the above result. 

The condition for the formation of a small po- 
laron roughly corresponds to Ep > t, where t is 
the nearest-neighbor transfer integral. According 
to the result Eq. (4), taking t = 50meV as repre- 
sentative for organic crystals and assuming a typi- 
cal value z <~ a (a = 7.2A in rubrene), we see that 
the formation of small polarons is expected as soon 
as P > 0.03. The physical parameters characteriz- 
ing the different dielectrics used in Rcfs. [2,3] are 
reported in Table I. 



3. Activated mobility 

The mobility of small polarons in the hopping 
regime can be evaluated by mapping the prob- 
lem onto a two-site cluster. [91 In this reduced 
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Table 1 

Physical properties of the gate dielectrics, ujq represents 
the typical range of phonon frequencies found in the litera- 
ture for the bulk materials, /3 is the dimensionless coupling 
constant. The dielectric constant of rubrene is ft = 3. The 
values reported for Si3N4 take into account the possible 
formation of the oxide SiO^N^ during device fabrication. 

model, the electrons interact with the long-range 
polarization field created by the interface optical 
phonons via an appropriate collective coordinate 



Q = E g (l - cos qa)M q 



i_ q ). In the adia- 



batic regime, where the lattice motion is much 
slower than the electrons, the electronic degrees of 
freedom can be integrated out leading to a double- 
well potential landscape for the coordinate Q. The 
mobility in the hopping regime can be expressed 
using standard techniques as 



ea 2 u> s 
h 27rT f 



-A/T 



(5) 



where the Kramers rate is determined by the fre- 
quency lo s of oscillation around the minima of the 
potential- well, and the activation energy A cor- 
responds to the height of the barrier. It is related 
to the polaron energy Ep defined in Eq. (4) by 



A = -fE P - t. 



(6) 



Here 7 is a parameter which depends on the shape 
of the interaction matrix element: [6] 



1 The prefactor in Eq. (5) is affected by several microscopic 
parameters such as the lattice geometry, the interaction 
with multiple phonon modes, the polaron size, dynamical 
friction, which can sensibly modify this result. 

2 In the above derivation we have implicitly assumed the 
adiabatic regime, which is valid when t 2 / (u! s ^/EpT) 3> 
1. A formula similar to Eq. (5) hods true in the oppo- 
site anti-adiabatic limit. In that case the activation bar- 
rier is A = fEp, and the prefactor must be replaced by 
(v^/3)(t 2 /VT5A). 
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Fig. 1. (reproduced from Ref. [3]) Left panel: polaronic part 
of the mobility vs. temperature in devices with different 
gate dielectrics. Full lines are fits according to Eq. (5). 
Right panel: extracted values of the activation barrier A 
vs. the dimensionless coupling /3. 
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(7) 



The second equality has been obtained using Eq. 
(2), which is valid for z > a. Note that for a local 
interaction (Holstein model, M q = const), from 
Eq. (7) one recovers 7 = 1/2. A calculation of 7 
including the effects of the discrete lattice is pre- 
sented in the Appendix. 

The temperature dependence of the mobility 
has been measured in Ref. [3] in devices using the 
different gate dielectrics listed in Table I. In prac- 
tice, the measured fi(T) includes contributions 
from other scattering mechanisms, such as the 
coupling with molecular vibrations inside the or- 
ganic crystal. These are entirely responsible for 
the mobility i^r(T) observed in devices in which 
the gate insulator is vacuum, and can be sub- 
tracted out from the raw data via the Matthiessen 
rule = fip 1 + /ip 1 . 

In Fig. 1 we report the polaronic part of the 
mobility, obtained with the above described pro- 
cedure. Upon increasing the dielectric constant of 
the gate insulator (which amounts to increasing 0) , 
the temperature dependence evolves towards an 
insulating-like behavior that can be fitted reason- 
ably well with Eq. (5). The values of the activation 
barrier A extracted from the fits are reported in 
the right panel. Despite the rather limited exper- 
imental temperature range, A is a linear function 
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of (3 for the devices with the strongest polarizabil- 
ities (S1O2, Al 2 03, Ta20s), as expected from Eqs. 
(4) and (6). Furthermore, from the slope of A vs 
(3, a typical distance z ~ 3A between the carriers 
and the interface can be estimated, which is com- 
parable with the lateral size of a rubrene molecule. 
Note that the effect of 7 < 1/2 [cf. Eq. (7)] was 
not considered in Ref . [3] , which led to a higher es- 
timate for this distance (z ~ 6.4A). 



4. Electron-electron interactions 

Having established that the electron-phonon in- 
teractions can be tuned by changing the polariz- 
ability of the gate material, we now investigate the 
possibility of controlling the Coulomb interactions 
between electrons. It is clear that these are negligi- 
ble at low carrier concentrations, i.e. in the linear 
regime of the Isd vs. Vq curves, where the conduc- 
ticity is directly proportional to the carrier density 
(see e.g. Fig. 1 in Ref. [3]). 

However, especially in the devices with strongly 
polarizable gate dielectrics, a sizable carrier den- 
sity can in principle be induced at the highest at- 
tainable gate voltages. For example, in a device 
with a Ta205 layer of thickness d = 500nm, a con- 
centration x ~ 0.15 holes per rubrene molecule can 
be reached at Vq = 100V. At such density, the ef- 
fects of the electrostatic repulsion between the car- 
riers cannot be neglected a priori. Here we propose 
a simple theoretical framework to describe their 
possible consequences on the device performances. 

Let us start by observing that the Coulomb in- 
teraction between two electrons in proximity to 
a polarizable dielectric is screened by the image 
charges as 



<p(r) = — 

K 



1 



1 



r yjr 2 + 4z 2 e s + n 



(8) 



In the following we shall assume that the average 
interelectron distance is larger than the distance z 
between the electrons and the interface, so that Eq. 
(8) can be replaced by its asymptotic expression 
2/( K + e s )/r. 

In a system of N interacting electrons, each 
charge feels instantaneously the sum of the repul- 



sive potentials of the others. Dynamical screening 
is unlikely in the regime under study, where the 
motion of the charge carriers is hindered by the 
formation of small polarons. For the same reason, 
if we focus on a given particle while it hops to its 
neighboring site, the relaxation of the remaining 
electrons can be neglected to a first approxima- 
tion: the collective rearrangement of N — 1 charges 
necessarily implies several hopping processes, and 
therefore occurs on a much longer time scale than 
the individual hopping event under consideration. 

From the above arguments, if the interacting 
fluid is initially at (or sufficiently close to) equilib- 
rium, it can be expected that moving a given par- 
ticle from a site to its nearest neighbor will have 
a net energy cost W, that adds to the polaronic 
hopping barrier A. In fact, such an energy cost can 
be easily included in the two-site model presented 
in the preceding Section. It leads to an increased 
value of the activation barrier A — > A+W/2, which 
eventually reduces the carrier mobility compared 
to the result Eq. (5) for independent polarons. PH 
The mobility in the presence of electron-electron 
interactions can thus be expressed as: 



KT) = MT)e- w ' 2T . 



(9) 



To get an order of magnitude estimate of this ef- 
fect, we replace the actual (uniform) distribution of 
polarons by a triangular array of point charges at 
the same density. This approximation is expected 
to be qualitatively correct for the following rea- 
sons: first of all, it is known[10] that the local cor- 
relations in an interacting charged liquid are very 
similar to the ones of a crystallized state. On the 
other hand, due to the long-range nature of the in- 
teraction potential Eq. (8) , the result for W should 
depend only weakly on the details of the charge 
distribution. Finally, the proposed lattice approx- 
imation enforces in a simple way the absence of 
relaxation of the electron fluid mentioned above. 

The potential energy of a given particle in a lat- 
tice can be expanded for small displacements u 
around its equilibrium position as: 



3 This expression is valid for W A. More generally, the 
activation barrier in a biased double- well potential is given 
by (A + iy/4) 2 /A. 
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AE lattlce {u) = (——J-^u 2 . (10) 

Here R s is the radius of the Wigner-Seitz disk 
n = (7ri? s )~ 1 , which is proportional to the average 
electron-electron distance, and£ = 0.8 is a geomet- 
rical factor characterizing the triangular lattice. In 
the present approximation, the additional energy 
cost for hopping between molecules induced by the 
electrostatic repulsion between the carriers is sim- 
ply given by W = AE lattice (a). Defining the car- 
rier concentration according to a/ R s = ^/tFx leads 
to 

2 e 2 

W = ( — (ttx) 3/2 = 0.3lx 3/2 eV (11) 

k + e s 2a 

in the case of Ta20s. For x — 0.15, we obtain a 
value W/2 ~ lOmeV that is consistently smaller 
than the polaronic barrier A ~ bbmeV . According 
to Eq. (9) the Coulomb interactions between the 
carriers can reduce the mobility by 30% at room 
temperature, and even more at lower temperature. 
Such behavior should be clearly visible as a bend- 
ing down in the Isd vs Vq curves at the highest 
attainable values of the gate voltage. [11,12] 

Note that the concentration that can be reached 
at a given voltage is proportional to the dielectric 
constant e s . As a result, for devices of equal thick- 
ness, the maximum value of W increases with e s 
as ~ e^J 2 /(k + e s ). The observation of Coulomb 
interaction effects is therefore more likely in FETs 
with highly polarizable gate materials. 



5. Concluding remarks 

In this paper we have reviewed a model which 
describes the interaction of the charge carri- 
ers with the interface phonons in single crystal 
organic FETs. Our analysis demonstrates that 
the electron-phonon interaction can be tuned by 
changing the dielectric polarizability of the gate 
material. This model consistently explains the 

4 It can be argued that such high values of Vq can modify 
the charge distribution along z, therefore modifying the 
value of the activation barrier according to Eqs. (4), (6) and 
(7). [5] This effect, however, should not play a dominant 
role if the polarons are already "small" . 



evolution of the carrier mobility from a metallic- 
like to an insulating-like behavior that has been 
recently observed in Refs. [2,3]: increasing the 
dielectric polarizability of the gate insulator re- 
sults in a crossover from the weak to the strong 
coupling regime where the carriers form small po- 
larons, which gives rise to a thermally activated 
mobility. Building on this model, the possibility 
of revealing experimentally the effects of the long 
range Coulomb repulsion between the carriers at 
high concentrations has been analyzed. 

In conclusion, although the main interest in or- 
ganic FETs comes from their potential applica- 
tions in " plastic electronics" , the opportunity of 
tuning several parameters such as the carrier den- 
sity, the electron-phonon and possibly the electron- 
electron interactions, makes them an ideal play- 
ground for fundamental physics. 



Appendix A. On the lattice cut-off of the 
electron-phonon interaction. 

Any linear electron-phonon interaction of the 
form of Eq. (1) can be transformed to real space as 

where g tj = J2 g M q e~ iqi - R ^ R i\ R t being the lat- 
tice sites. A simple prescription to define a discrete 
version of the Frohlich model, appropriate on a lat- 
tice, is to introduce a short-range cut-off of the or- 
der of the lattice spacing, which defines a new po- 
tential gf™* . The correct periodicity of the matrix 
element on the Brillouin zone can then be restored 
by setting 

M q = J2 eiqRi 9tf- (A.2) 
Ri 

It is clear that the new matrix element obeys 
M q+ c = M q with G any reciprocal lattice vector, 
which is not the case for M q of Eq. (2). By taking 
the continuous limit and letting the cutoff Rq — > 
one recovers M q — M q . 

Frohlich interaction in 3D. The essence of the 
Frohlich model is the interaction of a charge 
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(monopole) with a field of dipoles in an ionic ma- 
terial (+ and - ions). Correspondingly, the usual 
polar electron-phonon interaction in 3 dimensions, 
Mg D — Mo/q, scales in real space as 1/R 2 , namely 



M 



(A.3) 



The simplest way to introduce a lattice cut-off is to 
flatten off the potential at short range by setting 
9u = 9i,i+i- Alternatively, one can define 



3D, cut 



M 



2^{Rl+Rl) 



(A.4) 



with Ro of the order of the lattice spacing a. Trans- 
forming back g\f' cut to reciprocal space leads to 
the familiar Yukawa form. 

A physically sound choice for the cutoff length 
Ro in Eq. (A.4) is to set Ep equal to the energy 
of self-interaction of a uniform charge distributed 
on a cube of side a, [7] which yields Ro — 0.4a. 
However, the polarizability of the medium at short 
distances is certainly less than what is predicted 
by continuum theory so that the above estimate 
should be considered as a lower bound for Rq. 

Note that the potential Eq. (A.4) looks essen- 
tially similar to the one proposed in Ref. [8] to de- 
scribe the interaction with the apical oxygens in 
the superconducting cuprates, except for the dif- 
ferent decay at long distances, reflecting a differ- 
ent physical origin of the interaction (in that case, 
gij <~ 1/Rfj). As for the short distance cut-off, the 
authors of Ref. [8] use Ro — a. 

Polar interaction at interfaces. The interaction 
term of Eq. (2), that we rewrite here as M q = 
M e~ qz I \fq, already has a short distance cut-off 
determined by the distance z from the interface. 
Indeed, transforming to real space yields the fol- 
lowing Hypergeometric function 



9ij = 



M 



4 v / 7TZ 3 / : 



2^(3/4,5/4,1,-^./^), (A.5) 



that scales as l/R^j 2 for large Rij z, but tends 
to a constant gu — l/4y / 7rz 3 / 2 for R^ <C z. Again, 
a lattice cut-off can be introduced in this expres- 
sion either by flattening off the local term gu or by 
replacing z — > max(z, a). 



As a final remark, to avoid dealing with the spe- 
cial function 2-F1, a simpler model can be intro- 
duced that reproduces both its asymptotic limits 
(but slightly underestimates the polaron energy). 
It reads 



cMo 



9ij 



(A.6) 



where c = l/[4r(l/4)r(5/4)] ~ 0.076 and R = 
(40Fc) 2 / 3 z ~ 0.66z. 

The polaron energy in this model can be evalu- 
ated in real space as Ep = gfj/ui s . Performing 
the discrete sum on a square lattice in the case z = 
a gives a result which is about 30% smaller than 
the result of the continuous approximation. On the 
other hand, the value of the parameter 7 that de- 
termines the activation barrier increases by almost 
a factor 2, leading to an overall enhancement of 
the barrier. Although these numbers are quite sen- 
sitive to the choice of the cut-off prescription, they 
confirm that the continuous approximation used in 
the text is reasonably accurate for z > a. 
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